Experimental implications of quantum phase fluctuations in layered superconductors 
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I study the effect of quantum and thermal phase fluctuations on the in-plane and c-axis super- 
fluid stiffness of layered d-wave superconductors. First, 1 show that quantum phase fluctuations in 
the superconductor can be damped in the presence of external screening of Coulomb interactions, 
and suggest an experiment to test the importance of these fluctuations, by placing a metal in close 
proximity to the superconductor to induce such screening. Second, I show that a combination of 
quantum phase fluctuations and the linear temperature dependence of the in-plane superfluid stiff- 
ness leads to a linear temperature dependence of the c-axis penetration depth, below a temperature 
scale determined by the magnitude of in-plane dissipation. 

PACS numbers: 74.20.-z, 74.20.De, 74.72.-h, 74.25.Nf 



I. INTRODUCTION 



The linear temperature dependence of the low tem- 
perature in-plane penetration depth, X a b(T) is well es- 
tablished now in nearly all the high-T c superconductors 
jl], p. This behavior is most simply explained in terms 
of nodal quasiparticle excitations in a d-wave supercon- 
ductor ||. These low energy nodal quasiparticles then 
determine the low temperature behavior of thermody- 
namic properties and in-plane response functions. How- 
ever, since the high-T c superconductors (SC's) have a low 
superfluid density and a short coherence length compared 
to conventional superconductors, it is plausible that in 
addition to quasiparticle excitations, quantum and ther- 
mal phase fluctuations could also be important. The im- 
portance of such quantum phase fluctuations, especially 
for c-axis properties, has been emphasized earlier, from 
an analysis of the c-axis optical conductivity and sum 
rules d). 

Based on our study of an effective quantum XY phase 
action in the superconducting state, we had argued [|[ ||] 
that quantum phase fluctuations are important in the 
cuprate SC's at low temperature. With increasing tem- 
perature, in the presence of dissipation, the system grad- 
ually crosses over from a quantum fluctuation regime 
to a regime of classical phase fluctuations at a 

crossover scale, T c \ < T c . (This has been observed in 
recent experiments on thin films of conventional dirty s- 
wave superconductors]^] where T c \ ~ 0.94 T c .) It has 
been suggested ||, that this crossover scale T c \ <C T c 
in the high-T c systems, and thus classical phase fluctu- 
ations should dominate the low temperature supercon- 
ducting state properties. However, for parameter values 
relevant to the cuprate SC's, and even overestimating the 
magnitude of dissipation, which should lead to a smaller 
crossover temperature, we foundQ that T c \ ~ 20K. 
Since a linear T behavior of the superfluid stiffness has 
been observed in many high-T c SC's from ~ 3QK down to 
~ IK, with no evidence of such a thermal crossover, we 



conclude that the crossover temperature is high > 30K, 
and classical effects are irrelevant for in-plane proper- 
ties at low temperatures T <C T c . The low tempera- 
ture regime is thus described entirely by quasiparticle 
excitations around a ground state with strong quantum 
phase fluctuations. Recently, an experiment has been 
proposed JT(| to observe these quantum phase fluctu- 
ations by studying the excess current along the c-axis, 
after suppressing Josephson tunneling using an in-plane 
magnetic field. 

The presence of strong quantum phase fluctuations in 
the low temperature superconducting (SC) state raises 
several questions: Can the magnitude of low temper- 
ature quantum phase fluctuations be directly probed? 
What is the effect of such phase fluctuations on the c- 
axis superfluid stiffness? We shall study these questions 
in this paper. The principal results of this paper are as 
follows. (1) External screening of the Coulomb interac- 
tions in a superconductor damps quantum phase fluctu- 
ations. Based on this I suggest an experiment where a 
metal is placed in close proximity to the superconduc- 
tor to induce such external screening, similar to experi- 
ments which have been carried on gated Josephson junc- 
tion arrays pj. In the quantum fluctuation regime, rele- 
vant to our case, this screening reduces the fluctuations 
and leads to observable changes in the superfluid stiff- 
ness and its temperature dependence. In the classical 
regime, the screening would not lead to any change in 
the superfluid stiffness, since the fluctuations are already 
overdamped. Thus, this experiment would serve as a way 
to directly measure the magnitude of quantum fluctua- 
tions in the superconductor. (2) I show that quantum 
phase fluctuations can lead to a linear T increase of the 
c-axis penetration depth. This linear T slope arises from 
the linear T dependence of the in-plane superfluid stiff- 
ness and the Josephson coupling between layers. The 
magnitude of this effect is sensitive to the in-plane dis- 
sipation. The linear temperature dependence may be 
consistent with some low temperature c-axis penetration 
depth data pi pL §1 in YBa 2 Cu 3 7 _ 5 (YBCO) and 
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Bi2Sr 2 CaCu208+,5 (BSCCO) systems, thus serving as in- 
direct evidence of quantum phase fluctuations in the SC 
state. More generically however, the c-axis penetration 
depth does not appear to have a linear temperature de- 
pendence. I discuss possible implications of this towards 
the end. 



II. EFFECTIVE PHASE ACTION 

I begin by briefly reviewing the formalism of Refs 
U to study phase fluctuations in the low temperature 
superconducting state. We model the dynamics of the 
phase variables 0r(t) defined on a coarse-grained lattice 
(with lattice spacing equal to the coherence length, £o) 
by a quantum XY action. We have derived the following 
coarse-grained action for the layered d-wave SCs: 



■;r\vn\%XQ) 
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where d c is the interplane separation, £o is the in-plane 
coherence length, a = x,y,z and 7 y (Q) = (4 — 2 cos Q x — 
2 cos Q y ). The internal dissipation in the superconduc- 
tor arises from the conductivity of the electronic de- 
grees of freedom which have been integrated out, and 
is parametrized by a where, for simplicity, I assume 
a constant (ohmic) layer conductance ad c = (e 2 /h)a. 
The couplings J a are given by J x . y = J y = D®d c and 

J z = J ± = D° ± d c (£ /d c ) 2 , where D°,D° ± denote the 
bare in-plane and c-axis stiffnesses respectively, which 
are related to the penetration depths through A~ 2 ± = 

ATie 2 D r ^/h 2 c 2 . V(Q) = V(Q l{ a/£ ,Q ± ) denotes the 
scaled Coulomb interaction which arises from the coarse 
graining procedure, with 



V(Q) 



2ire 2 ad c 
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sinh(Q||d c /a) 



cosh((5||d c /a) — cosQj 



(2) 



being the Coulomb interaction appropriate for layered 
systems. Here ej, is some background dielectric constant, 
a is the original in-plane lattice spacing, and Qy, Q± are 
the in-plane and c-axis components of the momentum, 
measured in units of 1/a and 1 jd c respectively. 

I analyze the quantum XY action within the self con- 
sistent harmonic approximation jn| (SCHA). The SCHA 
replaces the above action by a trial harmonic theory with 
the renormalized stiffness D» I± chosen to minimize the 
free energy. It thus takes into account the renormaliza- 
tion of the stiffness due to anharmonic longitudinal phase 
fluctuations but ignores vortices and vortex-antivortex 
pairs (transverse fluctuations). The validity of this ap- 
proximation at low temperature, in the classical case 
where only n = Matsubara frequency is retained, has 
been confirmed in recent Monte Carlo simulations of the 



classical XY model ||. Carrying out the SCHA for our 
anisotropic case leads to 



D, 



D° exp(-« x )/2) 



(3) 



where S0 n , ± = (Ot,t — #r+a,r) with a — x,y for the in- 
plane (||) case and a — z for the c-axis (J.) case. The 
expectation values are evaluated in the renormalized har- 
monic theory, which yields 

K> - 2T f (0E7„(Q)G(Q,^) 
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where 7,,(Q) = (4 — 2cosQ x — 2cosQ a ) and 7 X (Q) = 
(2 — 2cosQ z ). Here, G(Q,w„) is the phase propagator 
given by 
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In what follows, I analyze the above equations in various 
limits, at low temperature, and also compare with a full 
self-consistent numerical solution Jl7[ of Eqs. (||) - (||). 



III. AN EXPERIMENTAL PROBE OF 
IN-PLANE QUANTUM PHASE FLUCTUATIONS 

In this Section, I study the damping of in-plane quan- 
tum phase fluctuations due to external screening of 
Coulomb interactions. This leads us naturally to a sug- 
gestion for an experiment to estimate the magnitude of 
these fluctuations in the cuprate SCs. Previous work 
on the damping of quantum phase fluctuations include 
studies on resistively shunted Josephson junctions arrays 
with short range charging energy jlq] and on a related 
model for superconductors with a low superfluid den- 
sity [ fl9| |. These works focussed on the quantum phase 
transition from an insulator to a superconductor, driven 
by increasing the strength of dissipation. Subsequently, 
such a quantum phase transition was observed experi- 
mentally in gated Josephson junction arrays 



If , where 



the quantum phase fluctuations were damped by exter- 
nal screening from a (gate-tuned) metallic bath. I show 
below within our phase fluctuation action, which differs 
in important respects from the earlier models, that such 
damping of quantum fluctuations through an external 
metallic bath leads to observable consequences for the 
low temperature superfluid stiffness. This may be used 
to study the importance of quantum fluctuations in the 
high-T c superconductors, with a set-up similar to the one 
used in the above experiment. 

Let us confine ourselves to the case of a thin su- 
perconducting film capacitively coupled to a metal- 
lic bath, and work in the two-dimensional (ID) limit 
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of the phase action. For electrons situated at an 
interface between vacuum, with a dielectric constant 
of unity, and a (metallic) substrate with e su b(w) = 
(too + 47ricr cx t/w), the Coulomb interaction is given by 
V2d(Q) = 47re 2 arf c / (1 + e su b(w)) Q» . Thus, the metal- 
lic bath provides dynamical screening of the Coulomb 
interactions [^9[ ^0), and leads to an additional term in 
the phase action in Eq. (Q), such that 

(a/27r) 7|| (Q)K| -> [a 7|| (Q) + a s Q„] (6) 

where a E = (fo/d c ) (<T cx td c ) / (e 2 / h) , and e h -> (l+£oo)/2 
in the Coulomb interaction in Eq. (||). Thus, external 
dissipation due the metallic bath (cr cx t) appears together 
with the internal dissipation from the electronic degrees 
of freedom of the superconductor which have been inte- 
grated out (a), and they can both lead to damping of 
quantum fluctuations in a similar manner. Note that 
our approach differs from the earlier work of Ref. Jl9| in 
two important ways, (i) We retain the dynamical term 
io 2 £ 2 d c /V(Q) in the action in Eq. ([!]), and the renormal- 
ized propagator in Eq. (||). (ii) We make a distinction 
between the dissipation which arises from the degrees of 
freedom internal to the superconductor and that from ex- 
ternal screening. This is reflected in the fact that the con- 
ductivities a and We appear with different Q-dependent 
coefficients in the above action pj| . 

To quantitatively study the effect of the external 
screening, let us consider parameter values relevant to the 
cuprate SC's. I choose Coo ~ 10, £o/a « 10 and W ~ 10 
as representative of the YBCO system p3| and any typ- 
ical substrate material. We model YBCO as a system of 
strongly coupled bilayers, with the bilayer stiffness be- 
ing twice the single layer stiffness, and an inter-bilayer 
spacing d c /a sa 3 being twice the mean layer spacing. 
This will simplify our calculation, since we do not have 
to introduce an additional parameter to distinguish intra- 
and inter-bilayer couplings; a more sophisticated calcu- 
lation would not lead to any qualitative or significant 
quantitative change. For the present analysis, we will 
use a bare bilayer stiffness and its linear T slope such 
that the renormalized stiffness J. (T) , for the above pa- 
rameters and with We = 0, leads to a penetration depth 
A t| (0) ~ 1600A and dXJdT w 4A/K in agreement with 
experiment 0|. We can then vary the external dissi- 
pation <7g, and study its effect on J,.(T). In Fig. |l|, I 
plot the behavior of the the stiffness J y (0), and its slope 
<2J|i /dT, for various values of We corresponding to differ- 
ing levels of external dissipation. 

It is clear that the stiffness J y (0) increases with increas- 
ing dissipation, as quantum fluctuations are damped, 
leading to a more ordered state. With increasing dis- 
sipation, (classical) thermal phase fluctuations also con- 
tribute to the normal fluid density, in addition to the 
'bare' quasiparticle contribution, leading to an enhance- 
ment in the slope of J,, (T) . A measurement of the stiff- 
ness and its temperature dependence in the presence of 




FIG. 1: The T = in-plane superfluid stiffness for various val- 
ues of external dissipation We- The inset shows the behavior 
of the linear T slope of (T). The above results have been ob- 
tained by solving the Eqs. (^)-(^) for a two dimensional case, 
assuming bare bilayer stiffness J» (0) ~ 150meV with a slope 
~ 0.65meV/K , and an internal dissipation W = 10. These val- 
ues have been chosen such that, for We = we obtain a pen- 
etration depth A (0) w 1600A, and a slope dX /dT « 
in agreement with experiment^. 



external screening would thus serve as a test of quan- 
tum phase fluctuations. This could be probed in low fre- 
quency optical experiments, or penetration depth experi- 
ments. A minor caveat concerning penetration depth ex- 
periments is that the effective (Pearl) penetration depth 
A c ff = X 2 /S where S is the film thickness, and hence one 
has cannot use the strictly 2D limit in the above calcula- 
tion but S has to be kept finite. This is not expected to 
lead to any qualitative changes in our results for J.. (T). 



IV. QUANTUM PHASE FLUCTUATIONS AND 
THE C-AXIS PENETRATION DEPTH 

Having shown that quantum phase fluctuations may 
affect the in-plane superfluid stiffness, I next turn to ex- 
amine the effect of these fluctuations on the c-axis su- 
perfluid stiffness. Many earlier studies of the c-axis stiff- 
ness have focusscd on the quasiparticle contribution to 
the normal fluid density, emphasizing the role of tunnel- 
ing matrix elements |23|] , or effects of disorder and pair 
tunneling on the temperature dependence of the critical 
current pj, |25|. Here, I study the effect of phase fluctua- 
tions on the c-axis superfluid stiffness within the SCHA, 
and find that phase fluctuations can lead to a linear tem- 
perature dependence of the c-axis superfluid stiffness and 
penetration depth. This linear T slope arises from the 
linear T dependence of the in-plane stiffness and the 
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Josephson coupling between layers, and it is a generic 
effect in such models. 

We shall begin by assuming that the bare c-axis stiff- 
ness is T-independent at low temperature. This as- 
sumption may be justified as follows. Quasiparticles 
may be important for in-plane properties, but the inter- 
plane tunneling matrix element p6[ being proportional to 
~ (cos k x — cos k y ) 2 , leads to a very small depletion of the 
superfluid density, since matrix element vanishes for in- 
plane nodal quasiparticles. In our calculation we shall 
neglect c-axis dissipation since the sharp Josephson plas- 
mon seen in experiments on BSCCO implies a very small 
c-axis conductivity, and the conductivity has been mea- 
sured to be small over a wide frequency range [ p7[ . In 
YBCO, there are added complications due to the chains, 
but I will ignore this. 

To study the effect of phase fluctuations, let us analyze 
the fluctuation integral for (SO 2 ^ ) in Eq. (EJ) for the physi- 
cally relevant case of J ± /J,, -C 1. Since we are interested 
only in J ± (T) we fix the in-plane stiffness, J,,(T), from 
experiment, and study the c-axis fluctuations (50 2 ± )(T). 
I will first present results for the case with no in-plane 
dissipation (a = 0) and later consider the effect of dis- 
sipation to see how these results are affected by a finite 
(7 =/= 0. Finally, I will compare our results with experi- 
mental data on BSCCO and YBCO. 



A. Dissipationless case (a — 0) 

For a = 0, we can first do the Matsubara summation 
in the fluctuation integral in Eq. (El). This leads to 



prefactor. Thus we are led to 



2 / t^t7 



V 



_„ (2nf ^V(J ||7|l +^7x) 



x cothj ^\/V(J n j n +J ± i ± 
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f2ne 2 
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where C\ is a constant of order unity, which depends only 
on (,o/d c , and which can be determined numerically for 
a given system. It is now easy to see that the linear T 
dependence of J, directly leads to a linear T dependence 
of (SO 2 ) and hence of J ± ,X ± . To relate the slope of X ± 
to the slope of J.. , we set J.. (T) = J (0) — aT. Using 
Eq. (0) we then get 



dX ± 
~dT 



= d 




(9) 



B. Non-zero dissipation (a 7^ 0) 

In the presence of dissipation, we can again analyze 
the Matsubara summation in the fluctuation integral to 
obtain the asymptotic low temperature behavior, similar 
to our earlier analysis for the in-plane fluctuations || . In 
this case, for a 3> 1 and T — > 0, we find 



{sol) 




27re 2 



1 



^ J\\ (T) 

d 3 Q 7_l7h 
_ x (2tt)3 ( 7|| J u + 7i Jj 2 

While the leading linear T dependence of (SO 2 ) 



(10) 



arises 



from the linear T dependence of J„ in the first term in 
the above equation, a T 2 dependence arises, primarily 
from the second term for a > 1. With J.. (T) = J.. (0) — 
aT, as before, the leading temperature dependence of the 
fluctuations is now given by 



where we have suppressed the Q dependence of -f ± | ,. and 

V. We next note that the temperature dependence of 
(SO 2 ) arises from two sources: (i) in the coth factor, 
which corresponds to thermal excitations of the plasmon 
mode, and (ii) in the prefactor, through the tempera- 
ture dependence of J.. . I have numerically checked that 
the coth factor may be set to unity at low temperature 
for the cases of interest, since the in-plane plasmon en- 
ergy is very large. Even when the c-axis Josephson plas- 
mon energy is very low, as in case of BSCCO, the phase 
space for this low energy excitation is small in the fluc- 
tuation integral and the characteristic plasmon energy is 
very high (see Fig. 1 of Ref. || ) . The crossover associated 
with this factor then only leads to large power laws for 
temperatures above the c-axis plasmon scale. Further, 
J ± /</,, C 1 means we can now safely set J ± = in the 



4T 



)(T)-(S6i)(0) = - 



2a 



where 



C 2 



d 3 Q 



a JA0) +C2 3 



7,i Tl 



T 



-* (2^) 3 (7 II +(JJ0)/J || (0)) 7 J^ 



(11) 



(12) 



is a constant depending on the anisotropy J ± / J., at T = 
0, which may be easily determined numerically. This 
leads to a crossover scale T 2 = 4a J y (0)/(C 2 (J 2 ) beyond 
which temperature the linear temperature dependence 
crosses over to a T 2 behavior. For T <C T2, we find a 
linear T behavior in X ± (T) with a slope 



dX ^ ' ! 



(13) 
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However, since the linear T to quadratic T 2 crossover 
temperature depends very sensitively on the dissipation, 
large dissipation might lead to a T 2 behavior down to 
the lowest observed temperatures. We therefore make 
estimates of this temperature scale X2 for BSCCO and 
YBCO6.95; we then compare our results for the temper- 
ature dependence of X x (T) with some experimental data 
on BSCCO and YBC0 6 . 95 . 



C. Comparison with experiments 

I begin by fixing the parameters connected to the in- 
plane stiffness in the phase action, and A^ (0) since we are 
only interested in our predictions for the temperature de- 
pendence of X ± (T) . I then compare the slope dX ± (T) / dT 
with existing experiments in BSCCO and YBCO sys- 
tems. 

For BSCCO, I use £ /a ~ 10 and d c /a ps 4, and set 
the bilayer stiffness J., (0) ps 80meV with its linear T 
slope a ps 0.8meV/ K . These values correspond [^8| to a 
penetration depth A y = 2100A with dXJdT ps IOA/K. 
I then set A^ (0) ps 150/im and choose a reasonable 
value (ct = 20) for the internal dissipation. This gives 
Ci « 0.6, C 2 ps 1.35, and T 2 ps 60K below which we 
expect to see linear T behavior; the low temperature 
slope dX ± /dT ps 0.075/im/ K is somewhat smaller than 
some experimentally reported values (l^, [ll| of ~ 0.25- 
0.3/xm/A'. 

For YBC0 6 .95, I use £ /a ps 10, d c /a ss 3.2, and 
set the bilayer stiffness J., (0) ps lOOmeF with its slope 
a ps 0.5meV/K. This leads to a penetration depth 
A y ss 1600A with dXJdT ps aX/K in agreement with 
experiment |l], Finally, I set A i (0) ps 1.1/xto, and use 
a reasonable value a — 10 for the intrinsic dissipation. 
This gives C x ps 0.5, C 2 ps 0.4, and T 2 > T c , below which 
we expect to see linear T behavior. Thus, the linear T 
behavior from phase fluctuations is expected to persist 
over a larger temperature scale in YBCO6.95. I then find 
the slope dX ± /dT ps 5A/K. This is somewhat smaller 



than some reports Jl3| of dX ± /dT ~ 15-20 A/if on this 
system. While this behavior was attributed jL3) to the 
effect of the chains in YBCO, the role of phase fluctua- 
tions could clearly also be important. The fluctuations 
at T = lead to a ~ 30% renormalization of the stiffness 
J ± (0), in reasonable agreement with an earlier c-axis con- 
ductivity sum rule analysis j| carried out for YBC06.6- 
More generically, the experimentally observed c-axis 
penetration depth in BSCCO and YBCO6.95 is reported 
to have a weaker temperature dependence (see Ref. JlJ] 
for instance), possibly ~ T 2 at low temperature. One 
possible reason for this discrepancy between the predic- 
tion of the quantum phase fluctuation model and the 
experiments could be the effects of disorder between the 
planes which deserves a more careful investigation. 

V. CONCLUSIONS 

We have seen that it is possible to directly probe the 
importance of in-plane quantum phase fluctuations in 
the high-T c superconductors for T <C T c through mea- 
surements of the in-plane penetration depth Am (T) in the 
presence of external screening. I have also shown that an 
indirect measure of quantum phase fluctuations may be 
obtained from studying the temperature dependence of 
the c-axis penetration depth, X ± (T). It is possible that 
disorder between the planes could affect our prediction, 
leading to a weaker temperature dependence more con- 
sistent with the experimental data; we leave this issue 
for future work. However the linear T behavior might 
still be observable in some clean materials. Experimen- 
tal tests of these would lead to a better understanding 
of the superconducting ground state and its low energy 
excitations. 
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